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ABSTRACT 

The identification of on-going planet formation requires the finest angular resolutions and deepest sensitiv¬ 
ities in observations inspired by state-of-the-art numerical simulations. Hydrodynamic simulations of planet- 
disk interactions predict the formation of circumplanetary disks (CPDs) around accreting planetary cores. 

These CPDs have eluded unequivocal detection -their identification requires predictions in CPD tracers. In 
this work, we aim to assess the observability of embedded CPDs with ALMA as features imprinted in the 
gas kinematics. We use 3D Smooth Particle Hydrodynamic (SPH) simulations of CPDs around 1 and 5 Mj^p 
planets at large stellocentric radii, in locally isothermal and adiabatic disks. The simulations are then connected 
with 3D radiative transfer for predictions in CO isotopologues. Observability is assessed by corrupting with 
realistic long baseline phase noise extracted from the recent HL Tau ALMA data. We find that the presence of 
a CPD produces distinct signposts: 1) compact emission separated in velocity from the overall circumstellar 
disk’s Keplerian pattern, 2) a strong impact on the velocity pattern when the Doppler shifted line emission 
sweeps across the CPD location, and 3) a local increase in the velocity dispersion. We test our predictions 
with a simulation tailored for HD 100546 -which has a reported protoplanet candidate. We find that the CPDs 
are detectable in all 3 signposts with ALMA Cycle 3 capabilities for both 1 and 5 Mj^p protoplanets, when 
embedded in an isothermal disk. 

Subject headings: protoplanetary disks 


1. INTRODUCTION 


Planets are expected to form during the evolution o f cir¬ 
cumstellar disks of gas and dust (e.g., Armitage||2011 ). As 
a protoplanetary core grows to approximately a Saturn mass. 


it becomes massive enough to open a gap in the disk ( |Lin 
|& Papaloizou||1986t Lubow et al.|1999| ). This process is the 


result of competition between gravitational, viscous and pres 
sure torque s exerted onto the disk by the planet and by the 
disk itself (|Crida et al.||2006|). The gap splits the disk into 
two radially distinct zones. There are multiple detections of 
dust-d epleted gaps and cavit i es in well studied p rotoplanetary 
disks ( [Andrews et al.|[20TT| [Perez et ^|2015| ), but interest¬ 
ingly, no unambiguous detection of a forming planet has yet 
been confirmed. 

Hydrodynamical models of planet-disk interactions show 
that a single giant protoplanet continues accreting fro m the 
outer disk at formidable rates (^2xl0~^ Mj^pyr" ^, [Gres 


sel ( 


et al.||2013[[Shabram & Boley|2013j[Szulagyi et al.jf2(JlT|) 

even after its gap is evacuated ( [Kley|1999[[Papaloizou & Nel-[ 

son[ [2005[ ). Accretion streams converge onto the vicinity of 


the giant developing a circumplanetary disk (CPD) through 
which ang ular momentum disperses thus regulating planetary 
accretion ( [Lubow et al.|1999 [Ayliffe & Bate|2009] ). 

In two-dimensional (2D) simulations, strong shocks appear 
near the planet’s Hill sphere leading to excessive infiow redi- 
rect ed towards the planet, rapidly depleting t he CPD mate¬ 
rial ( [Lubow et al.|1999t[D’Angelo et al.|2002 ). These shocks 
appear much weaker in three-dimensional (3D) calculations, 
leading to more persi stent CPD structures ( Bate et aL][2003 


[Ayliffe & Bate[[2009[ ). Grid-based simulations show that the 


inflow of gas onto the protoplanet mostly happens in the ver¬ 
tical direction, a llowing material to cross the shock front nea r 


2010[ [Szulagyi et al 
VIHD simulations by 


ressel 


the Hill radius ( [Machida et al. 

Similar results were found for 
[etar[ ( |MT3] ). 

Most models show th at CPD radii truncate to about 1/3 
labram & Boley[[20T3]) or even 1/2 


(Ayliffe & Bate[[2QQ9[ [Shabram & Bol^ 

( Gressel et al. |2013{ ' Szulagyi et al.|2014[ ) 


of the planet’s Hill 
radius. This implies that a Jupiter mass planet on a 100 AU or¬ 
bit could bear a CPD with a diameter of 4.5-7 AU. If located 
at 100 pc away from Earth, the projected diameter translates 
into ^45-70 mas, within range of modern astronomical in¬ 
strumentation. 

Gas-giant protoplanet candidates have been detected em- 


bedded in the HD 100546 and HD 169142 disks, two in each 

system (Quanz et al.|20131 

[Reggiani et al.|2014[ Currie et al. 

2014 Biller et al.|i014|i. B 

loth systems are young and nearby 


circumstellar disks with confirmed gaps ( [Osorio et al.[[2014 
Walsh et ^[2014[ ). The two directly imaged detections are 


sperez @ das .uchile. cl 


L' (3.8/im) emission blobs: one is at 0.5” angular separation 
(^52 AU at 100 pc) from HD 100546 (also detected in M' 
emission at 4.8/im, [Quanz et al.[[2014| , and the other is at 
0.16” (^23 AU at 145 pc) from HD 169142. The compact 
source in HD 169142 falls w ithin a symmetric ga p imaged 
in polarized scattered light by [Quanz et al. (2013b). Interest¬ 
ingly, there are no near-IR counterparts for either candidate, 
supporting the idea that these are in fact accreti ng gas gian ts 
with SEDs driven by circumplanetary accretion ( Zhu|2015 ). 

These CPDs have eluded unambiguous detection , mainly 
bec ause o f the lack of predictions on CPD tracers. [Isella et 
al. ( 2014[ ) presented deep continuum observations to detect a 
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Table 1 

Simulations parameters 


Parameter 

Value 

Rin 

0.35 Rp 

Rout 

1.85 Rp 

Initial S profile 

S oc R-i(*) 

Initial H/R profile 

Circumstellar disk masses^ 

0.05 (constant with Ri) 

Isothermal (SPHl, SPH2) 

5 X lO-^M© 

Adiabatic (SPH3) 

7 X W-^Mq 


Note. — (*) normalized so that Rp = 1 AU would give E = 100 g cm ^ at 
R = Rp. (t) for the nominal disk model presented in Sec.l^ 


CPD around a protoplanet candidate in LkCa 15, with no pos¬ 
itive results. CPD continuum emission may in fact be scant 
since large dust grains (> 100 /im), probed by sub-mm contin¬ 
uum observations, are sieved out of the planet-induced gap by 
the outer disk pressure bump. Interestingly, only small grains 
(<100/im) enter the dust-depleted cavity ( |Zhu et aL||2012| ) 
and make it into the CPD. This dust-filtration is evidenced by 
the nu merous dust-depleted cavities seen in sub-mm observa¬ 
tions ( [Andrews et al.|20lT| . CPD observational signposts to 
unambiguously confirm forming planet candidates are scant. 
Realistic 3D simulations coupled with radiative transfer are 
needed to predict the observational signatures of circumplan- 
etary disks at various wavelengths. 

In this work, we connect 3D Smooth Particle Hydrodynam¬ 
ics (SPH) of CPD dynamics (Section with 3D radiative 
transfer (Section |3]) in the context of interferometric obser¬ 
vations of line emission. We aim to study the observability of 


CPDs as features imprinted in the gas kinematics (Section H). 
We compute predictions tailored for the Atacama Large Mi - 
limeter Array (ALMA) for common gas tracers (Section [43 ). 
Implications and conclusions are discussed in Section 


2. NUMERICAL SIMULATIONS 

We carried out a set of 3D SPH simulations to characterise 
the CPD morphology and kinematics, which we feed into a ra¬ 
diative transfer code (see Section [^. We address the question 
of observability signposts of a single snap CPD embedded in 
its parent disk, after the simulations have reached a reasonably 
steady state. Exploring how planet-disk interactions evolve 
with time is beyond the scope of this paper. 


2 . 1 . Three-dimensional SPH Simulations 

We used a mod ified version of the SPH code Gadget- 
2 ( |Springel||2Q05| ) to perform a suite of 3D simulations of a 
planet embedded in a protoplanetary disk. The code has been 
modified to make it more suitable for sim ulating planets em¬ 
bedded in disks (see jDunhill et al.|2013| ). We have simulated 
the full disk azimuth but a restricted radial range. 

We implement radia l bound ary conditions in a manner sim¬ 
ilar to I Ay liffe & B ate j ( [20 1 Qj ) , reducing the range of the sim¬ 
ulated aisF3EIIe3tIIl providing high resolution around the 
planet. Although spiral waves launched by the planet refiect 
off the boundary, they do not affect the disk dynamics at the 
radius of the planet. A fuller description of this method will 
be given in a forthcoming paper (Dunhill et al., in prep). We 
summarise used disk parameters in Table 1. 

To reduce the effect of transient waves on the planet at the 
start of the simulation, we include an initial gap in the disk us¬ 
ing the prescription of ( jLubow & D’ Angelo|2006| ). We model 
the planet as a point mass potential with a sink radius, inside 


which any gas particles are swallowed and their mass and mo¬ 
mentum added to the planet. The potential is unsoftened out¬ 
side the sink radius. Initially, i^sink = 0.05 and decays 
exponentially to i^sink = 0.001 after approximately 8 or¬ 
bits. This corresponds to ^2 Jupiter radii for Rp = 1 AU. 

We performed three different runs, two locally isothermal 
(SPHl and SPH2, where T{R) is enforced so that H/R re¬ 
mains constant) and one with an adiabatic equation of state 
(SPH3, using an adiabatic index 7 = 5/3). It has been shown 
before that using an isothermal EOS yields slightly larger cir- 
cumplanetary disks th an using more realistic radiation physics 
([Ayliffe & Bate|2009|). SPHl has a planet mass Mp = 1 Mj, 
while SPH2 ana SPH3 have Mp = 5 Mj, all orbiting a 1 Mq 
star. 

Eor the simulations with Mp = 5 Mj^p, we model the disk 
using 2 million SPH particles. In order to achieve an equiva¬ 
lent resolution within the CPD, we used 16 million SPH par¬ 
ticles in SPHl. We vertically resolve the CPD into ^ 4 SPH 
smoothing lengths h (typically h ^ 10“^i?p in the CPD mid¬ 
plane) e nsuring that w e do not under-estimate the midplane 
density ( |Nelson|2006| ). At this resolution, the artificial viscos¬ 
ity in the simulations gives an effective Shakura-Sunyaev al¬ 
pha parameter of a ^ 0.005 within the CPDs. Due to the ex¬ 
treme computational expense, we halted SPHl after 10 orbits 
of the circumstellar disk, although lower-resolution tests indi¬ 
cate that it has settled into a steady state by this time. SPH2 
and SPH3 were halted after 50 orbits (see Fig. 1). 

These simulations are limited in that we neglect complex ra¬ 
diation physics, including only viscous and shock heating but 
not passive heating. However, the CPD structures we focus 
on are still present in rad iation hydrodynamic runs as shown 
by jAyliffe & Bate| ( [2009 1. 

3. RADIATIVE TRANSFER PREDICTIONS 


The main driver of this investigation is to study under which 
conditions an accreting protoplanet would be detectable 
through ALMA observations of line emission. We have cho¬ 
sen bright and commonly observed CO transitions which lie 
within the sub-millimeter (sub-mm) range in ALMA. Rota¬ 
tional transitions of and the isotopologues, and 

C^^O, are k nown to trace the gas in protoplane tary disk gaps 
and cavities ( Bruderer|2Q13t[Perez et al.|2Q15 ). Most impor¬ 
tantly, these lines contain important kinematic information, 
essential to detect companion objects embedded in the gas in¬ 
side dust-depleted cavities. 

We consider two model disks for our analysis. A nominal 
disk located at 100 pc, inclined by 20 ° and hosting a planet at 
100 AU, used to illustrate the observational features revealed 
through CO kinematics. The second model is tailored for the 
HD 100546 disk, with an inclinati on of 42° and a 5 Mj^p 
planet candidate at Rp = 52 AU (jQuanz et al.||2014|). Our 
simulations are resampled accordingly m Cartesian c oordi¬ 
nates using a linear interpolation scheme via SPLASH ([Price 
2007|). The Cartesian cells are perfectly cubic, with a cell size 


of 0.014 Rp, which, after scaling becomes 1.4 AU and 0.7 AU 
for Rp = 100 AU and 52 AU, respectively. 

For the isothermal runs (SPHl and SPH2), temperature is 
an imposed function of orbital radius, while for the adiabatic 
run (SPH3) we use the SPH internal energies to calculate tem¬ 
perature by assuming a standard mean molecular weight and 
adiabatic index. We scale the discs using a T{R) oc 
consistent with measurements of flari ng discs ( [Kenyon & 


Hartmann[[1987[ Andrews et al.[[MTT] ). After scaling, tem- 
peratures at 100 AU reach ^60 K for the isothermal disks and 
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Figure 1. Surface density maps after 10 and 50 orbits for the isothermal and adiabatic simulations listed in Tabled Axes are in planet orbital radius Rp units. 
Left and middle show isothermal runs for 1 and 5 Mj^p planet candidates. Right. Adiabatic simulations for a 5 (50 orbits). Upper right insets show a 

close up to the CP D kinematics. The planet bears a disk and drives a spiral wave without disrupting the disk structure heavily. Figures made using the SPLASH 
code jPrice|2007) . 


-2000 K for SPH3, well above CO freeze-out (20-25 K). 

In the event of CO freeze-out, a similar RT calculation can 
be applied to species with enhanced abundances where CO is 
depleted, such as DCO+ or N 2 H+, or species that are formed 
by surface reactions with CO ice, such as H 2 CO, which have 
recently been detected in disks at or bey ond the CO snow-line 
( |Qi et al.|2013a|bt [Mathews et al.|2013t . 

We compute synthetic images in CO(2-l) with the radiative 
transfer code RADMC3D (Dullemond et al. 2015). Line ra¬ 
diative transfer is done in LTE, using molecular data from the 
LAMDA databasfl Line widths include thermal broadening 
and a local (spati^y unresolved) microscopic turbulence set 
to a constant value of 0.1 kms“^. We used a fiducial molec¬ 
ular abundance relative of H 2 relative to ^^CO of 1( )~^. We 
adopted an ISM i sotopic abundance ^^C /^^C of 76 (|Stahl et 


al.|2008| ) and 500 ( [Wilson & Rood] 1994| ) for ^^CO andC^^O, 


respectively. 

Channel maps are rendered using RADMC3d ray-tracing. 
The results are synthetic data cubes centered on the star, with 
a total width in velocity of —16 km and individual chan¬ 
nels of 0.1 km s“^. These data cubes represent our sky model 
which is subsequently Fourier transformed and resampled to 
ALMA’s visibility plane. 

For completeness, we calculate the continuum assuming a 
simple dust di stribution consisting of 30 per-cent amorphous 
carbon grains ( |Li & GreenbergH 19971 ) and 70 per-cent astro¬ 
nomical silicates~ fDraine & Lee[l984j i, following the gas den¬ 
sity. Grain size distribution follow a power-law with expo¬ 
nent 3.5 from 0.05 to 1000 jim. We compute dust opacities 
with Mie theory. CPD dust continuum predictions are ad¬ 
dressed elsewhere in th e literature (see. Wolf & D’Angeloj 
2005||Isella et al.|2014|), but are strongly a ffected by gas pres¬ 
sure bumps ( D’Angelo & Podolak||2015| ). Continuum emis¬ 
sion in sub-rnm observations arises froni dust thermal radia¬ 
tion and does not contain kinematic information, thereby, it is 
not within the scope of this paper. 


4. RESULTS AND DISCUSSION 

Fig. [^presents our ^^CO channel maps predictions, at ve¬ 
locities -1.0, 0.0, -1-1.0, -1-2.0 and -1-3.0 kms“^, for the nominal 
disk model. Top and middle panels are isothermal simulations 
for 1 and 5 Mjup planets (SPHl and SPH2, respectively), 
while bottom panels show adiabatic results (SPH3). 


^ http://www.strw.leidenuniv.nl/-moldata/ 


We find that the presence of a CPD produces deviations 
from circumstellar Keplerian kinematics. The circumplane- 
tary disk is in itself a miniature Keplerian accretion disk em¬ 
bedded in the gap, and it can be separated in velocity from 
the overall Keplerian pattern of the circumstellar disk. This 
velocity separation can be seen in Fig. as compact emission 
at the CPD position. This compact emission is persistent over 
a velocity range given by the CPD kinematics. 

Spectra extracted from the vicinity of the CPDs produces 
a broader profile when compared with a spectrum extracted 
from the point symmetric opposite side of the disk (see 
Fig.|^. The profile extracted from the vicinity of the 1 Mj^p 
planet shows a broad core on top of even broader wings 
(top spectrum in Fig. [^. The middle panel in Fig. shows 
that the isothermal CPD around the 5 Mj^p planet exhibits 
a distinct double-peaked profile, with a peak separation of 
— 1.5 kms~^. Full CPD spectra span over Au —2.5 kms~^ 
for the isothermal 1 Mj^p CPD, while the 5 Mj^p covers 
more than /S.v —4 kms~^. The double-peaked profile for 
the 1 Mjup case is unresolved at 0.1 kms“^ resolution. The 
width of these spectral features may inform on the size of their 
respective CPDs, and ultimately on the mass of the accreting 
planet via their Hill radii. Fig. also shows that the CPD line 
wings for SPHl and SPH2 ena in an abrupt shoulder. This 
is likely due to our planet accretion model; a point mass with 
a sink radius inside which particles are swallowed and their 
kinematics cannot be sampled. 

The adiabatic CPD spectrum does not reveal distinctive fea¬ 
tures (bottom spectra in Fig.[^. The compact CPD emission 
is also much less distinctive in the adiabatic disk (bottom pan¬ 
els in Fig. 1^. In opposition to the isothermal disk where all 
compressive work is radiated away immediately, the adiabatic 
disk cannot cool, reaching temperatues of —2000 K around 
the CPD. This causes the hot gas to rapidly fill back the gap 
preventing a CPD to fully develop (see second moment map 
in Fig. 1^, which hinders clear identification of spectral fea¬ 
tures in an adiabatic fiow. Indeed, adiabatic and isothermal 
disks represent two extrema of the phenomenon that we are 
modelling. 

A second signpost of planet formation arises when the 
Doppler shifted line emission of the circumstellar disk’s Ke¬ 
plerian pattern sweeps across the CPD location. The but- 
terfiy pattern becomes strongly bent and twisted, while the 
point-symmetric location at the opposite side of the disk re¬ 
mains undisturbed. This can be seen in the central panels of 
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Figure 2. Predictions for emission based on SPHl, SPH2 and SPH3 (top, middle and bottom, respectively). The crosses show the central star and circles 

indicate the CPD position. From left to right, the maps velocities correspond to -1.0, 0.0, +1.0, +2.0 and +3.0 km . Channel widths are 0.1 km . Maps at 
0.0 km represent systemic velocity. Inset frames show predictions for ^^CO. Color scale and contours are logarithmic. See Sectionj^for a full description. 
Fluxes are given in Jy pixel“^, where each synthetic pixel is 12 mas. 


Fig.J^ The insets show the same twisted pattern but for ^^CO. 
Optically thicker than ^^CO, the maps still reveal the 

kinematic bend, even for a 1 Mjup planet whose gap is shal¬ 
lower and its appear optically thick around the CPD. 

For SPH3 the CPD vicinity is much hotter than in the isother¬ 
mal cases, producing enhanced ^^CO emission (see bottom 
inset in Fig.[^. 

As noted m Section 3, our choice of temperature profile 
in SPHl and SPH2 is inconsistent with the scaling used for 
the RT calculations. Self-consistency should result in a CPD 
approximately twice as thick as in the SPH, reducing the mid¬ 
plane densities and emitted fiuxes by a similar factor. It is pos¬ 
sible that the signatures highlighted in Figures 2 and 3 would 
thus be at lower contrast to the background, but still present 
at detectable levels as the velocities are largely unaffected. 
CO freeze out at the CPD’s location would be prevented by 
adding a backgro und temperature to acco unt for accretion ra¬ 
diation feedback ( [Montesinos et al.|2Q15|) and incident ra dia- 
tion from their environment (see Shabram & Boley|2013| ). In¬ 
cluding self-consistent thermal physics in future simulations 
will settle this discrepancy. 


4.1. HD 100546 through a 15 km baseline sub-mm 
observation 

To assess the observability of these kinematic CPD sign¬ 
posts, we performed a second calculation tailored for the pro¬ 
toplanet candidate in HD 100546, based on the SPH2 run. We 
tied the fiuxes of our model to match approximately previous 
C O obse rvations of this source ( |Walsh et al.|2014| [Pineda et] 
aL]|2014 ). We filtered our sky model using the i4i;-coverage 


from the long base line (^15 km) Science Verification obser¬ 
vations of HL Tau ( [Partnership et al. |20 15 j ) . We corrupted our 
model with thermal and phase noises extracted directly from 
the HL Tau dataset. The simulated observation was then self- 
calibrated and CLEANed using routines in the CASA pack¬ 
age. 

Fig. ffl shows the HD 100546 simulated ALMA observa¬ 
tion. Lot and right panels show selected channels for CO 
and ^^CO emission. The upper panels illustrate the recov¬ 
ery of the bent Keplerian locus, while lower panels show 
how the CPD compact emission can also be detected. Both 
emission features are recovered at the 5 a level (rms noise is 
1 mJy beam ^). 

5. CONCLUSIONS 


The presence of a CPD produces distinct signposts in simu¬ 
lated CO maps: a striking compact emission separated in ve¬ 
locity from the overall Keplerian pattern of the circumstellar 
disk; a strong infiuence on the velocity pattern of the gas when 
the Doppler shifted line emission sweeps across the CPD lo¬ 
cation; and a local increase in the velocity dispersion. More¬ 
over, for the locally isothermal simulation with a 5 Mjup, the 
CPD spectra even exhibits a double-peaked profile. These dis¬ 
tinctive features rely on kinematics and can reveal the pres¬ 
ence of an embedded CPD perturber even in optically thick 
tracers like ^^CO. 

The feasibility of an ALMA observation of HD 100546 was 
assessed by corrupting our synthetic visibilities with realis¬ 
tic phase and thermal noises extracted from the HL Tau long 
baseline campaign. We found that these CPDs are detectable 
in all 3 signposts with ALMA Cycle 3 capabilities for both, 1 
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Figure 3. CPD spectra extracted from an 80 mas aperture (in radius) 

centered on the planet (circle in Fig. [^. Top, middle and bottom spectra 
correspond to SPHl, SPH2 and SPH3, respectively. Squared data points (grey 
curve) are spectra extracted from the point symmetric location at the opposite 
side of the disk. Insets show second moment maps (velocity dispersion, a) 
calculated over the region of interest. 


and 5 Mj^p protoplanets, when embedded in a locally isother¬ 
mal disk. On the other hand, in the pessimistic case of an 
adiabatic disk the CPD formation was hampered, hence de¬ 
tectability is scant. Previous radiation hydrodynamic calcula¬ 
tions, show that CPDs are better described by isothermal disks 
and that the adiabatic simulati on is in fact a rather pesimistic 
case (see |Ayliffe & Bate[2009 their fig. 12). 

The immediate vicinity of the planet’s Hill sphere, in¬ 
cluding the CPD itself, offers an environment for gas-phase 
physics which produces distinctive kinematic observational 
features. Future ALMA long baseline observations of gas 
tracers could detect these signposts of planet formation, and 
provide not only confirmation of forming planets but also 
valuable kinematic information on CPD physics. 
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Figure 4. Channel map signpost predictions for HD 100546’s CPD candi¬ 
date. Calculation based on the 5 Mj^p isothermal SPH2 run for two chan¬ 
nels centered at systemic velocity (bottom panels) and +2 km/s (top pan¬ 
els). Continuum subtracted ^^CO and ^^CO(2-l) emission after corrupting 
by HL Tau’s phase and thermal noise, a) and b) show the twisted Keplerian 
butterfly pattern, while c) and d) illustrate the striking CPD emission offset 
from the global Keplerian pattern. Channel widths are binned to 0.5 km s“ ^. 
The rms noise is ~1 mJy beam“^ in each bin 

57, and QUIMAL-130008 provided funding for improve¬ 
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